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Abstract 

We perform a comprehensive analysis of target mass corrections (TMCs) to spin-averaged 
structure functions and asymmetries at next-to-leading order. Several different prescriptions for 
' TMCs are considered, including the operator product expansion, and various approximations to 

o," it) collinear factorization, and ^-scaling. We assess the impact of each of these on a number of ob- 

. servables, such as the neutron to proton F2 structure function ratio, and parity-violating electron 

scattering asymmetries for protons and deuterons which are sensitive to 7Z interference effects. 
The corrections from higher order radiative and nuclear effects on the parity-violating deuteron 
asymmetry are also quantified. 
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I. INTRODUCTION 



Tremendous progress has been made in recent years in the quest to reveal the structure of 
the nucleon at its deepest levels. Traditionally deep-inelastic scattering (DIS) of leptons has 
been the primary tool used to study nucleon structure at large values of the four-momentum 
transfer squared, Q^, where the nucleon's quark and gluon (or parton) constituents can be 
cleanly resolved. Here the theoretical tools are well developed, and the nucleon's structure 
can be conveniently parametrized in terms of universal longitudinal momentum distribution 
functions of individual quarks and gluons. Global analyses of DIS and other hard processes 
have been highly successful in correlating data over a wide range of kinematics, producing 
fits of parton distribution functions (PDFs) to next-to-leading order (in the strong coupling 
parameter as) accuracy and beyond 1-^. 

While the perturbative domain of large and small parton momentum fraction x has 
received considerable attention both experimentally and theoretically, the region of large x 
and low (~ 1 — 2 GeV^), where nonperturbative effects play a greater role, has been 
relatively poorly explored. This is perhaps not too surprising given the difficulty in reliably 
computing the various corrections that are needed to describe data in this region. Examples 
of nonperturbative effects that are relevant here include target mass corrections (TMCs) 
associated with finite values of M^/Q^, where M is the nucleon mass, higher twist terms 
arising from long-range nonperturbative multi-parton correlations, and nuclear corrections 
in experiments involving deuterium or heavier nuclei, which are important at large x for any 

The large-x region has been particularly difficult to access experimentally, especially in 
high-energy colliders, due to the rapidly falling cross sections as x — )■ 1. The most ex- 
tensive data set available that covers this region has been from experiments at SLAG j^. 
More recently, progress on this front has been made with DIS structure function mea- 
surements at Jefferson Lab, utilizing the high luminosities and duty factors available with 
the GEBAF accelerator. Indeed, an impressive body of very high-precision data has now 
been accumulated over the last decade on various structure functions, including accurate 
longitudinal-transverse separations needed for model-independent determinations of the F2 
and Fl structure functions of protons and nuclei [l^ . Future plans at the energy upgraded 
Jefferson Lab involve extending the DIS measurements to even larger x (x ~ 0.85) with 
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planned experiments Ill4l3l| to measure the ratio of to m quark distributions, as well as 



search for effects such as charge symmetry violation in PDFs and tests of the Standard 
Model in parity-violating DIS asymmetries. 

The new data have the potential to provide strong constraints on PDFs at large x, 
where currently uncertainties remain significant. Several recent analyses Is, 7| have in fact 
attempted to utilize data at low values of and invariant final state hadron masses = 
+ (5^(1 — x)/x. Stable fits of leading twist PDFs could be obtained 6| for W down to 
~ 3 GeV and ^ 1.7 GeV^, as long as TMC and higher twist corrections were accounted 
for. Aside from its intrinsic value, better knowledge of PDFs at large x may also be important 
for searches of new physics signals in collider experiments such as at the Tevatron or the LHC 



at large rapidities or for heavy mass particles 



14j |. as well as at more central rapidities where 



uncertainties in large- a; 



small X and large 15 1. 



'DFs at low can, through evolution, affect cross sections at 



The increased kinematic reach of the future high-precision DIS measurements calls for 
a careful evaluation of the relevant nonperturbative corrections in order to unambiguously 
extract information on leading twist PDFs or new physics signals. The effects that are most 
amenable to direct computation, in principle, are the TMCs. As discussed by Nachtmann 



161], these effects are in fact associated with leading twist operators (hence contain no addi- 



tional information on the nonperturbative parton correlations), even though they give rise 
to Q'^ jv^ = 4x^M^/Q^ corrections, where u = Q'^/2Mx is the energy transfer. Nachtmann 
further showed that one could generalize the standard operator product expansion (OPE) 
of structure function moments to finite such that only operators of a specific twist would 
appear at a given order in 1/Q^. The resulting target mass corrected structure functions can 
hen be derived through an inverse Mellin transformation, as shown by Georgi and Politzer 



17l | (for a review of TMCs in the OPE approach see Ref. 18|). 



Later an alternative formulation in terms of coUinear factorization (CF) was used by 

n 

Ellis, Furmanski and Petronzio [IQi] to derive TMCs including the effects of off-shell partons 
and parton transverse motion. While the OPE and CF formulations yield identical results 
for leading twist PDFs, they differ in the details of how the target mass corrections are 
manifested at finite Q^. Other versions of TMCs were subsequently derived 2Ch22|| within 
the CF formalism using various assumptions about the intrinsic properties of partons and 
higher twist contributions, leading to rather large differences in some cases [22']. Some 
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of the phenomenological implications of the different TMC prescriptions were discussed in 



Refs. 
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22| . including differences between leading order and next-to-leading order (NLO) 
results; however, the effects on observables have not been systematically investigated. We 
do so in this paper. 

In Sec. ini we summarize the main results for TMCs in the OPE and various CF formu- 
lations for the Fi, F2, F3 and Fl structure functions at NLO, and illustrate the differences 
numerically. Implications for various observables are discussed in Sec. Illlt including the 
ratio of neutron to proton F2 structure functions, which constrain the d/u PDF ratio at 
large x, longitudinal to transverse cross section ratios R, and parity-violating (PV) DIS 
asymmetries on the proton and deuteron which are sensitive to 7Z interference structure 
functions. We also quantify the effects of perturbative NLO corrections on the ratio 
for the 'jZ interference, about which nothing is known empirically, and of nuclear effects on 
the deuteron PV asymmetries. Some finite-Q^ effects on PV asymmetries were investigated 



previously in Ref. 23|, U^, and higher-twists in deuteron PV asymmetries in Refs. 251-129 1. 



Finally, in Sec. IIVI we draw some conclusions and outline possible extensions of this work. 

II. TARGET MASS CORRECTIONS 

In this section we review the kinematic corrections to structure functions arising from 
scattering at finite values of jv^ . We consider several frameworks for the TMCs, including 
the conventional one based on the operator product expansion, and various approximations 
to it, as well as a number of prescriptions using collinear factorization at leading and next- 
to-leading order in a^. The structure functions for the scattering of an unpolarized lepton 



from an unpolarized nucleon are defined in terms of the nucleon hadron tensor 



30|, 



= ^ / d'ze^-{v\{Jl{z).JM\\v) (la) 



H / \ H / \ / p ■ q 

- ie^ua^q p"— , lb 

2p ■ q 

where is the electromagnetic or weak current operator for a given virtual boson (7, Z or 
W^). Here p and q are the nucleon and exchanged boson four-momenta, respectively, with 
q' = -Q'. 
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The structure functions Fi^2 are related to the product of two vector or two axial-vector 
currents, while F3 arises from the interference of vector and axial-vector currents. The Fi 
structure function is proportional to the transverse virtual boson cross section, and F2 is 
given by a combination of transverse and longitudinal cross sections. It is convenient to also 
introduce the longitudinal structure function, 

Fl{x, Q^) = p'F^ix, Q') - 2xF,{x, Q') , (2) 

where 



P' = l + ^^- (3) 



Ax^M^ 

'Q 

In the following we will summarize target mass corrections for each of these structure func 
tions computed within the various approaches outlined above. 



A. Operator product expansion 



Target mass corrections to structure functions were first systematically considered by 



Georgi and Politzer 



171] in the framework of the operator product expansion. Here the 



twist-two quark bilinears in the product of currents J^^Ju in Eq. ffTal) are modified with the 
introduction of covariant derivatives, ip'^^D^^'^ ■ ■ ■ D'^^ip; since each derivative D^' increases 
both the dimension and spin of the operator by one unit, the twist (dimension minus spin) 
remains unchanged. The expansion in terms of covariant derivatives yields a series in 
with coefficients given by moments of structure functions. The resulting target mass cor- 
rected structure functions are then accessed through an inverse Mellin transformation, which 
gives 



17| 



^OPE, 



nOPE, 



4p 

+ 



2p2 



4p2 
3x(p2 - 1 
2p' 

X(p2 - 1) 



p' 



2xp 



p'-l 



p- 

jp' - 1: 

2p3 



2xp 



2xp 



(4a) 
(4b) 
(4c) 
(4d) 



where Fj^'^'' are the structure functions in the M'^/Q'^ — )■ limit, evaluated at the modified 



scaling variable ^ [16|, |31| 



2x 



(5) 



which approaches x as M'^/Q'^ — )■ 0. The functions /i2, 92 and are associated with higher 
order terms in M'^/Q'^ and are given by [3] 



du ■ 



du / (if- 



f,g2 



J du{u-^)- 



1 ^(0) 

du ^ 



Fr(u,Q') 



u 



(6a) 
(6b) 
(6c) 



(Note that the function 92 here should not be confused with the spin-dependent 92 structure 
function measured in polarized lepton-nucleon scattering.) 

The expressions in Eqs. (jl]) are known to suffer from the "threshold problem" , in which 
the target mass corrected (leading twist) structure functions do not vanish 1, and 

are in fact nonzero in the kinematically forbidden x > 1 region, where for a proton target 
baryon number conservation would be violated. This is clear from the 0{1) terms in Eqs. (jl]) 
in which the massless functions F^'^^ are evaluated at ^. Because at any finite value one 
has ^ < = C,{x = 1) < 1, for any input function f/'^^ which is nonzero for < a; < 1, 



the target mass corrected function at x = 1 will not vanish, 



OPE, 



X 



1,Q2 < 00) > 0. 



A number of attempts have been made to ameliorate the threshold problem 
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33| using 



various prescriptions and ansdtze, although none of these is unique and without additional 



complications 



18|. 



Recently, Kulagin and Petti 3J] showed that by expanding the target mass corrected 
structure functions to leading order in l/Q"^, the resulting functions have the correct x 1 
limits, 



\{5-p')Fi'\x,Q') - l(p^-l) [xFr'(x,g^ 



{x,Q' 



1 



(2 - p') (x, Q') - -Jp'- 1) [xF^'ix, 



- h2ix,Q^) 
6xh2{x, 



(x,g2 



1 



(p'-l) xFf (x,Q2)- 4xh2{x,Q^) 



i(7-3p^)Ff(x,g^) - i(p^-l) 



xFf (x,g^)- 2/i3(x,g' 



(7a) 
(7b) 
(7c) 

(7d) 



While avoiding the threshold problem, this prescription, however, raises the question of 
whether the 1/g^ approximation is sufficiently accurate for structure functions near x ~ 1 



at moderate . To test the convergence of the expansion at large x, we further expand 
the OPE results (j4]) to include 0{1/Q^) corrections, 



16 ^ ^ ^ 16a; ^ 



(0), 



16 ^ V >v ; -r 22 1 



^'^"'x,Q^) - -h,{x,Q') + ^g,{x,Q' 
4 8a; 



^a) 



23 
16 



F?\x,Q') + 



3x 



x,Q') 



+ ^F^'^"{x,Q') - 3xh,{x,Q') + l92{x,Q') 



(8b) 



Fl"'\x.Q') + (p^-1)^ 
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+ ^F, 



(8c) 



\X) 



+ 



F3^/«^(x,g^) + (p^-i)^ 



^Ff(x,g^)+ ?^Ff'(x,Q^) 



13 
16 



16' 



(8d) 



where the first (-F/°^') and second (Fj^^") derivatives of the structure functions are with 
respect to x. In fact, one can show that for a structure function that behaves at large x as 
(1 — x)", the target mass corrected result will vanish in the a; — )■ 1 limit up to order l/Q^""^ 
in the expansion. For n ^ 3, as is typical for nucleon structure functions, the threshold 
problem will therefore appear only at order 

The accuracy of the l/Q^ expansions is illustrated in Fig. [H where in order to isolate the 
target mass effect from the specific form of the structure function parametrization we have 
taken for simplicity the form F2 ~ (1 — x)^. Both the l/Q^ and l/Q^ approximations are 
found to reproduce the full OPE result very well up to x ~ 0.6, but significant deviations 
are visible at larger x. Furthermore, while there is a modest improvement in the agreement 
with the exact result for 0.6 ^ x ^ 0.8 after inclusion of the terms, both expansions 

appear to break down for x ^ 0.8. The reliability of a low order l/Q"^ expansion is therefore 
questionable at these x values, and hence their efficacy in removing the x — )■ 1 threshold 
problem. 

Since the integrals in the functions /;.2,3 and g2 can be time consuming to compute nu- 
merically, Schienbein et al. p^] found phenomenological analytic forms which approximate 
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FIG. 1: Ratio of the target mass corrected F2 structure functions using the 1/Q2 (sohd, red), 
(long-dashed, green) and phenomenological (short-dashed, blue) OPE approximations compared 
with the exact OPE result, Eq. ()4bp . Note that the phenomenological OPE approximation is 
almost indistinguishable from the exact OPE result, while the and expansions deviate 

from this for x ^ 0.6. 



the target mass corrected F2 and F3 structure functions in Eqs. ( l4b|) and fj4d|) by 



2p3 

(1+p) 



p(l + p) 

ip' - 1) , 



(9a) 
(9b) 



2p2 2p(l + p) 

These turn out to be rather good approximations to the exact results, as Fig. [1] illustrates 
for the F2 case. For all values of x, the phenomenological approximation ( l9a|) stays within 
5% of the full OPE resuh. 



B. Collinear factorization 



An alternative approach to TMCs relies on the collinear factorization (CF) formalism 



19l422|. which makes use of the factorization theorem to relate the hadronic tensor for 
lepton-hadron scattering to that for scattering from a parton. Here parton distributions 
are formulated directly in momentum space, avoiding the need to perform an inverse Mellin 
transform to obtain the PDF from its moments. An advantage of the CF formalism for 
TMCs is that it can be extended to other hard scattering processes, such as semi-inclusive 
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DIS [35], where an OPE is not available. 



1. Ellis, Furmanski and Petronzio 

The first study of TMCs within CF was made by Ellis, Furmanski, and Petronzio (EFP) 
191], who analyzed the virtual photon-hadron scattering amplitude using a Feynman diagram 
technique to expand the hard scattering term about the coUinear direction, incorporating 
both parton off-shellness (or interactions) and parton transverse momentum in twist-4 con- 



tributions 



36(1 . Using the same notation as for the OPE TMCs above, the EFP results for 



the target mass corrected structure functions are given by 

1 + p + 

Fr{x,Q') = (^'Q') + 7^^^2(e,Q^), (lOb) 

Fnx,Q') = Fi'\^,Q') + ^^^h,{^,Q^), (10c) 

Fn^,Q') = lFt\C,Q') + J^^h.i^^Q") > (lOd) 

where again the Pj^'^^ refer to the uncorrected structure functions, and /i2,3 are given in 
Eqs. (Note that the definition of the longitudinal structure function in EFP differs from 
the usual definition ([2]) by a factor x, and the F2 structure function is proportional to what 
EFP call the "transverse" structure function, which in standard usage is proportional to Fi.) 
Because the massless functions f/"'' are evaluated at ^, the target mass corrected structure 
functions will suffer from the same threshold problem as in the OPE analysis in Eqs. (jl]). 



While the expressions in Eqs. f IlOp were derived in Ref. 19| at leading order in ag, in this 
work we will assume their validity also at NLO. 

The prefactors for the leading terms proportional to Fj^^^ in Eqs. ffTOj) are remarkably 
close to those for the leading terms in the OPE expressions in Eqs. (j4]). To first order in 
1/Q^, the leading term prefactors for Fi in both OPE and EFP reduce to (1 — x'^M'^/Q'^). 
Similarly, the F2 prefactors both reduce to (1 — Ax"^ M"^ / Q"^) , while those for F^ reduce to 1. 
For the F^ structure function, however, the (9(1/(5^) prefactor is (1 — 3a;^M^/Q^) for OPE, 
whereas for the EFP CF resuh it is (1 - 2x^M^/Q'^). 

At leading order in the massless limit the longitudinal structure function vanishes iden- 
tically. At NLO, however, it receives contributions from both quark and gluon PDFs con- 



voluted with the respective hard coefficient functions. For electromagnetic scattering, for 
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example, one has 

where c' = 2 e^, and ^2^^^'^^^ is given by the leading order expression for F2^'^\ Similar 
expressions hold also for the longitudinal structure functions associated with other elec- 
troweak currents. In our numerical calculations discussed below we will always compute 
at NLO. 

It is important also to note that Eqs. ffTOj) have been derived considering Feynman dia- 
grams with 2 or 4 legs attached to the hadronic correlator (see Figs. 2 and 3 of Ref. 19|). 
which for M = give rise to twist-2 and twist-4 contributions to the structure functions, 
respectively. For M ^ 0, however, the quark and gluon equations of motion allow one to 
extract a twist-2 contribution from the 4-leg diagrams, which when added to the twist-2 tar- 
get mass correction yields the full result in Eqs. ffTOl) . It is an interesting question whether 
by resumming the twist-2 parts of n-leg diagrams one would be able to recover the TMC 
expressions ( fTOi) . 



2. Accardi and Qiu 

In both the EFP and OPE treatments of TMCs, the resulting structure functions are 
nonzero for x > 1. The analysis of Accardi and Qiu (AQ) [22] traced this problem to baryon 
number nonconservation in the handbag diagram for M ^ 0. Working with 2-leg diagrams 
only, in contrast to EFP who also consider 4-leg diagrams up to twist-4, the AQ target mass 
corrected structure functions are given by 22 1 



Fi^Q(x,Q2) = Fi'\^,Q'), (12a) 

i^2^^(x,g^) = (e,Q^), (12b) 

Ft^'M') = ^^°He,Q^), (12c) 

F3^Q(x,Q^) = -Ff (e,Q^). (12d) 



P 



Here the functions f/"'' are defined as 



/ 



z 
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where C/ are the perturbatively calculable hard coefficient functions for a given parton flavor 
/, including parton charge factors, (pj are the parton densities of the nucleon, and the sum 
is taken over all active flavors. The upper limit in Eq. f[T^ ensures that the target mass 
corrected structure functions vanish for x > 1, as required by kinematics, although jet mass 
corrections need to be introduced in order to render the target mass corrected functions 

.....„.H ,. .,.„......„„,....„.„„... 

twist-4 diagrams as in Ref. [19| in order to establish a more direct correspondence between 
the AQ and EFP approaches. Of course, for M^/Q^ — i- the upper limit of integration in 
Eq. (IT^ i s 1, and both approaches recover the standard factorization theorem for structure 



functions 



39|. 



3. ^-scaling 

When the upper limit of integration in Eq. ( 1131) is taken to be 1, the AQ structure 



functions reduce to the simple ^-scaling (^-S) form introduced by Aivazis et al. [20|] and 
used by Kretzer and Reno [21,]. The target mass corrected structure functions in this case 
are simply given by 

Fl-\x,Q^) = Fi'\^,Q'), (14a) 
Fl-'M') = ^Ff(e,Q^), (14b) 

Ffi^^Q') = ^-^FjP\^,Q-), (14c) 
Ft'{x,Q') = -^Ft\^,Q'). (14d) 

Note that the form of the target mass corrected functions in Eqs. ( JT4l) closely resembles that 
in Eqs. (fT2|) . with the two forms equivalent at leading order. At this order the structure 
functions satisfy a modified Callan-Gross relation {22 1, 

p'Ft\x,Q') = 2xFt\x,Q'). (15) 

The leading order ,^-scaling structure functions are also related to the leading, 0{1) terms 
of the OPE expressions in Eqs. (jl]), 

^OPE(lead,ng)^^^ Q^) = F^^, Q') , (16) 
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where the prefactor, to order l/Q"^., is given by (1 — x^M'^/Q'^). In fact, the ^-scahng 
formulas f|T^ would coincide with the EFP results in Sec. IIIB II in the absence of 4-leg 
Feynman diagrams [l9[ |. 



C. TMC comparisons 

The effects of the different TMC prescriptions on structure functions are illustrated in 
Figs. [2] - O for the F^, -F^, F]^ and structure functions of the proton, respectively. 

(The results for structure functions associated with other boson exchanges, such as W~ , Z, 
or 7Z interference, are very similar to these.) The uncorrected proton structure functions 



f/^^ are constructed from the CTEQ- Jefferson Lab (CJ) global PDF fits 
at = 2 GeV^. For each of the structure functions the effects of TMCs become more 
prominent with increasing x, and naturally their magnitude decreases at larger Q^. 



evaluated 



o 




w 
O 




X X 

FIG. 2: Ratio of target mass corrected to uncorrected (left) or to OPE (right) F'^ proton structure 
functions at = 2 GeV^ for the OPE (solid, red), EFP (short-dashed, blue), ^-S (long-dashed, 
green), and AQ (dot-dashed, orange) TMC prescriptions. Note that the OPE and EFP results are 
almost indistinguishable for x < 0.6. 



For the F^ structure function in Fig. O the deviation from unity of the ratio of target 
mass corrected to uncorrected functions ranges from ~ 10% at x = 0.4, to more than a 
actor 2 for x ^ 0.7. The model dependence of the TMCs to F^ is relatively weak; the OPE 



2Q|, 



2l| and AQ 



17 1 and EFP [19 [results are similar to within a few percent for all x, while the ^-scaling 



22| prescriptions differ from the OPE by ^ 5% and 15%, respectively, for 



X 



^ 0.8. (Results for x ^ 0.9 are not shown as the input nucleon PDFs are not constrained 
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in this region and display numerical instability at a; ^ 0.95.) In fact, at low and moderate x 
the OPE and EFP TMCs track each other rather closely, as expected from the equality of 
their leading term prefactors at order l/Q"^. Similarly, the AQ and ^-scaling prescriptions 
are much closer to each other than to the OPE and EFP results, as may be anticipated from 
the structure of the respective TMC expressions in Eqs. (IT^ and f[T^ . 




w 
o 




FIG. 3: As in Fig. [21 but for the F2 proton structure function. 

Qualitatively similar behavior is seen also for the target mass corrected F2 structure 
function in Fig. [31 Here a dip in the ratio of corrected to uncorrected functions at a; ~ 0.4, 
however, delays the sharp rise above unity to x ^ 0.6. As for F^, the EFP result agrees 
with the OPE to a few percent over the entire x range, and the AQ and ^-scaling ratios 
are almost identical for x < 0.4. The two sets of ratios differ by ^ 7% for x < 0.7, before 
diverging somewhat as x — i- 1. 

For the Fl structure function in Fig. [H the differences between the various TMC pre- 
scriptions are more dramatic. The OPE and EFP ratios begin to rise steeply at low x, with 
the corrected functions exceeding the uncorrected ones by more than a factor 5 already by 
X = 0.5. The AQ and ratios, on the other hand, rise above unity at much higher x, 
reminiscent of the Fi ratios in Fig. [2l The two sets of corrections differ by more than 50% 
by X = 0.3, and by ^ 80% for x > 0.8. This qualitatively rather different behavior can be 
understood by directly comparing Eqs. and ([TU]) to Eqs. ([T^ and (IT^ . Unlike the ^-S 
and AQ prescriptions, the OPE and EFP F^ results include terms involving integrals over 
F2, which is generally ^ Fl. In fact, to leading order with no TMCs, the Fl structure 
function vanishes, and adding NLO corrections within the AQ or ^-S prescriptions does not 
produce a significant increase. In contrast, the OPE and EFP prescriptions always receive 
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FIG. 4: As in Fig. [21 but for the proton structure function. 

large F2 contributions, making the target mass corrected to uncorrected ratio considerably 
larger in these approaches. 

The strong correlations between the OPE and EFP predictions are not as visible for 
the F3 structure function, which, unlike the other structure functions, differs already at 
0{1/Q'^). The general shape of the TMC ratio, illustrated in Fig. Elfor the F^"^ structure 
function, resembles that for in Fig. |3l but with a rise above unity beginning at lower x. 
The various prescriptions agree to ~ 10% for x ^ 0.4, and ~ 40% for x ^ 0.8, but generally 
display more spread than in F'^ or F^. 
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FIG. 5: As in Fig. [21 but for the -^3^^ proton structure function. 
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III. IMPLICATIONS FOR OBSERVABLES 



Having examined the differences between the various TMC prescriptions on individual 
structure functions, in this section we consider the effects of TMCs, and in particular their 
model dependence, on several observables that will be measured in upcoming experiments. 
These include the ratio of the neutron to proton F2 structure functions, the longitudinal to 
transverse (LT) cross section ratios, and parity-violating deep-inelastic scattering asymme- 
tries for the proton and deuteron. 



A. Neutron to proton ratio 

Historically, the ratio of (i to u quark distributions in the proton has been determined 
primarily through the ratio of neutron to proton electromagnetic DIS cross sections, 

d^a^P/dxdy FT \l-y- [p2 - 1 - 2pV(l + Wp)] /4 J ' ^ ^ 

where y = u/E is the fractional energy transfer from the lepton, and W'^ is the ratio of the 
longitudinal to transverse cross sections, or structure functions, 

= _L (18) 

for nucleon A^. With the assumption that = W^, the ratio of cross sections becomes 
the ratio of F2 structure functions, Rnp — )■ FT / FT ■ To leading order, the ratio f lT7|) is 
then given by Rnp = {1 + ^d/u) / {A + d/u), which illustrates the sensitivity to the d/u PDF 
ratio. In practice, differences between R'^p and R'^"' generated perturbatively at NLO have 
a negligible effect on the ratio Rnp at the kinematics considered here. 

The absence of free neutron targets has meant that in practice inclusive deuterium struc- 
ture function data has been used to obtain indirect information on the neutron, and hence 
the d quark. This procedure is known to suffer from significant model dependence at large 



values of x 



m 



4l|, leading to several novel new experiments being proposed ll|-|l3| to 
determine the d/u ratio with minimal nuclear model uncertainties. In order for these mea- 
surements to be unambiguously analyzed, it is important to quantify the extent of TMC 
uncertainties at the kinematics of the experiments, which will typically reach a maximum 
g2 ^ 10 - 15 GeV2 &tx^ 0.8. 
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FIG. 6: Ratio of target mass corrected [Rnp] to uncorrected (Rnp) neutron to proton F2 structure 
function ratios at = 2 GeV^ (left) and = 10 GeV^ (right), for the OPE (solid, red), EFP 
(short-dashed, blue), S^-S (long-dashed, green), and AQ (dot-dashed, orange) TMC prescriptions. 
The shaded band represents the uncertainty in the ratio d/u for the OPE result. Note that the 
AQ and ^-S results are almost identical. 



The effects of TMCs on the ratio Rnp are illustrated in Fig. [6] for the CJ PDFs 







at 



Q"^ = 2 GeV^ (left) and 10 GeV^ (right). The shaded bands represent the d/u uncertainty 
range as applied to the OPE TMC prescription, with the central solid (red) curve denoting 
the median value for OPE calculated with the same PDFs as the other TMC prescriptions. 
The target mass corrections at = 2 GeV^ are sizable, reaching ^ 25% — 30% at x = 0.7 
for the OPE and EFP prescriptions, and ~ 12% for the .^-S and AQ results. At the higher 
g2 = 10 GeV^ value the TMCs decrease to ^ 5% and ^ 2% for the OPE/EFP and ^-S/AQ 
calculations, respectively. Treating each of the TMC prescriptions on equal footing, this 
would suggest an uncertainty due to TMCs of ^ 3% for all values of x accessible in the 



planned experiments H 



13|. 



The TMC uncertainty can be compared with the range of Rnp predicted from PDFs 
extracted under different assumptions about the size of nuclear corrections in deuterium, 
which currently represents the largest uncertainty in the d/u ratio at a; ^ 0.5 [7]. This is 
illustrated in the bands in Fig. |6l which represent the Rnp ratio evaluated from the range of 
CJ PDFs [7] for the OPE TMC prescription. The results show that for = 2 GeV^ the 
uncertainty resulting from nuclear corrections is some 2 — 3 times larger than that associated 
with TMCs at X = 0.8. Both the uncertainties in the d quark PDF and in the TMCs decrease 
as X decreases, albeit more slowly for the latter. At x = 0.6, in fact, the two uncertainties 
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FIG. 7: Neutron to proton F2 structure function ratio iinp at = 2 GeV^ (left) and = 
10 GeV^ (right), with the shaded bands representing the d/u uncertainty for uncorrected structure 
functions (sohd, red) and the two extremal TMCs, EFP (fong-dashed, green) and AQ (short-dashed, 
blue). 

are comparable, while for x ^ 0.4 the TMC uncertainty is actually larger. 

With increasing both the TMCs and their uncertainties decrease, while the uncertainty 
in the leading twist PDFs remains approximately independent, as can be seen more 
clearly in the comparison of the absolute values of Rnp in Fig. [71 At = 2 GeV^ the 



bands representing the Rnp ratio evaluated from the CJ PDFs 7] using different TMC 
prescriptions (specifically the extremal EFP and AQ results) do not overlap until x ~ 0.75, 
meaning that at smaller x the true d/u behavior will be obscured by the relatively large 
TMC model uncertainty. Interestingly, the TMCs actually decrease the nuclear uncertainty 
range at lower Q^, since the action of the x C, rescaling is to effectively feed information 
from lower x in the uncorrected functions (which have relatively small d/u uncertainty) 
to higher x values (where the d/u uncertainty is larger). Consequently, at higher the 
sensitivity to the d/u ratio increases both due to the smaller spread of results for different 
TMC prescriptions, and to the weakening of the TMC effect in moving strength from lower 
X for a particular TMC prescription. This is indeed visible in Fig. [7] for the Rnp ratio at 
Q2 = 10 GeV^ in which the EFP and AQ extremal TMC bands very nearly coincide over 
the entire x range, as well as with the ratio computed without TMCs. Such values of 
will therefore be required in order to cleanly extricate the d/u ratio from measurements of 
the neutron to proton ratio without ambiguities associated with TMC model dependence. 
Finally, we also note that in global QCD fits of PDFs it was recently found [6] that the 
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impact of the model dependence of TMCs on leading twist PDFs is reduced significantly 
with the inclusion of a phenomenological higher twist term in the F2 structure function 
parametrization, with the two effects partially compensating each other. Since the higher 
twist contribution to is independent of that for F2, a similar cancellation may be expected 
also when including data in the global fits. The uncertainties in extracted PDFs induced 
by the model dependence of TMCs may therefore be smaller than those suggested in Fig. [6] 
if the data are analyzed within a global PDF context. 



B. Longitudinal to transverse structure function ratios 



While the LT cross section ratio W'^ is expected to play a relatively minor role in the 
measurements of the neutron to proton F2 structure function ratio in Eq. f|T7|) . mostly be- 
cause of the cancellation between the proton and neutron values, the effects of TMCs 
on the ratio itself may be more significant. This was already suggested by the large prescrip- 
tion dependence of TMCs for the longitudinal structure function Fl in Fig. HJ The effects 
of TMCs on the LT ratio are also important to quantify in connection with establishing the 
low-Q^ behavior of at finite x, to determine the onset of gauge invariance constraints 
on the longitudinal structure function 42 1. 





FIG. 8: Longitudinal to transverse cross section ratio for the proton at = 2 GeV^ (left) 
and Q2 = 10 GeV2 (right), for no TMCs (double-dashed, black), the OPE (soHd, red), EFP 
(short-dashed, blue), ^-S (long-dashed, green), and AQ (dot-dashed, orange) TMC prescriptions. 



In Fig. [H] we ihustrate the TMC effects on E?p for Q"^ = 2 and 10 GeV^ for each of the 
TMC prescriptions considered. All of the TMCs increase the magnitude of the ratio, 
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with the AQ and ^-S prescriptions having a relatively modest effect (approximately a factor 
2 for a; ^ 0.6 - 0.8 at = 2 GeV^, but only a few percent at = 10 GeV^), while the 
EFP and OPE both altering the ratio significantly for x ^ 0.1. The enhancement of for 
the latter is predicted to be about an order of magnitude for x ~ 0.6 — 0.8 at Q"^ = 2 GeV^, 
and still a factor of 3 - 4 at = 10 GeV^. 

Some differences are also expected between the longitudinal to transverse cross section 
ratios at NLO for processes involving electromagnetic and weak currents. In particular, as 
will be discussed in more detail in Sec. lIIICl below. asymmetries measured in parity- violating 
electron scattering are sensitive to interference effects between 7 and Z boson exchange, and 
differences between the W and LT ratios can affect the measured asymmetries 23l. |24|. 



OS 





FIG. 9: Ratio of K^^ to LT ratios for the proton computed at NLO for = 2 GeV^ (left) 
and Q2 = 10 GeV^ (right), for no TMCs (double-dashed, black), the OPE (solid, red), EFP 
(short-dashed, blue), ^-S (long-dashed, green), and AQ (dot-dashed, orange) TMC prescriptions. 



In Fig. inithe ratio of the proton to W LT ratios is shown at = 2 and 10 GeV^. 
While at leading order both of these ratios are zero, at NLO the different relative contri- 
butions from quark PDFs to the electromagnetic and 7Z interference structure functions 
leads to deviations of the ratios from unity of up to ~ 4% at = 2 GeV^, and up to ~ 2% 
at = 10 GeV^. The effects of the TMCs are again very small for the ^-scaling and AQ 
prescriptions, but more significant for the OPE and EFP results. Overall, the spread in the 
TMC predictions for the /W ratio amounts to ^ 4 — 5% for x between 0.6 and 0.8 at 

= 2 GeV^, and < 2% at Q"^ = 10 GeV^. Note that the dip in the ratios at a; < 0.1, which 
is insensitive to TMCs, reflects the greater role played by gluons at low x, but is mostly 



irrelevant for the kinematics of the proposed experiments 
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FIG. 10: As in Fig. H but for the deuteron to W LT ratio. 



For the case of the isoscalar deuteron target, stronger cancellations between the quark 
content of and W are expected to lead to smaller deviations of their ratio from unity 
at large x. This is indeed observed in Fig. (TDl where again the dip in the ratio at very 
low X is associated with NLO gluon dominance of the LT ratios as x — )■ 0. At x = 0.2, 
for example, the gluonic content of Fl suppresses the deuteron / W ratio by ^ 2% for 

= 2 GeV^, and ~ 1% for = 10 GeV^. At higher x the deviations decrease until the 
ratio approaches unity asymptotically as x — )■ 1. In the region of x where the LT ratios are 
dominated by quarks, the fact that the same isoscalar combination of quark PDFs enters 
both the electromagnetic and interference structure functions leads to almost negligible 
TMC effects. The absence of significant TMC effects in the deuteron ratio is, as expected, 
even more clearly visible at the higher value. 

Finally, for experiments involving deuteron targets one needs to account for the fact 
that the nucleons in the nucleus are bound and hence have structure functions that differ 
from those of free nucleons. In Fig. [TT] the LT ratio for the 7Z interference structure 
functions of the deuteron is shown relative to that for a free isoscalar nucleon (defined as 
proton + neutron) at = 2 GeV^, using the nuclear smearing model of Refs. |3J, |43| (see 
also Ref. [^). The shape of the deuteron to nucleon ratio computed at NLO in the absence 
of TMCs displays a dramatic rise above unity with increasing x that is characteristic of the 
nuclear EMC ratio [44]. The effects of Fermi motion in fact lead to a divergent ratio at 
X = 1. The inclusion of TMCs suppresses the rise at large x, from ^ 60% at x = 0.6 to 
20% for the ^-scaling and AQ prescriptions, and to ~ 5% for the OPE and EFP cases, with 
larger differences as x — )■ 1. This suppression arises because TMCs shift strength in Fl from 
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FIG. 11: Ratio of the NLO LT ratios for the deuteron and isoscalar nucleon {N = p + n) at 
Q2 = 2 GeV^ for no TMCs (double-dashed, black) and the OPE (sohd, red), EFP (short-dashed, 
blue), ^-S (long-dashed, green), and AQ (dot-dashed, orange) TMC prescriptions. 



small X to large x, thereby lessening the impact of the nuclear smearing. Since the TMC 
effects in the OPE and EFP prescriptions are larger than that for the ^-S and AQ cases (see 
Fig. H]), the nuclear corrections for the former in Fig. [11] are smaller. 

The smearing corrections to the electromagnetic LT ratio W are almost identical to those 
in Fig. [m Consequently, the net effect on the ratio W'^ / W for the deuteron computed with 
or without nuclear corrections is ^ 0.05%, and can be neglected for the kinematics of interest. 



C. Parity- violating DIS 

Measurements of parity-violating deep-inelastic scattering (PVDIS) asymmetries on the 
proton have been proposed at Jefferson Lab 13| to provide independent constraints on 
the ratio of d/u quark distributions at large x, free of the nuclear correction uncertainties 
associated with deuterium measurements [451]. In the case of deuteron targets, where much 
of the dependence on hadron structure effects cancels, PVDIS asymmetries are sensitive to 



several effects, .ndudrng charge symmetry violafou i,r PDFs M. o. to Standard Mode. 

parameters whose precise measurement can reveal signals for new physics [25|, |47[ . In this 
section we examine the effects of TMCs on the PVDIS asymmetries of the proton and 
deuteron, and discuss the phenomenological implications of their uncertainties on future 
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planned experiments. 

The PV asymmetry is defined through the difference and sum of the inclusive cross sec- 
tions for scattering either a right-handed (R) or left-handed (L) electron from an unpolarized 

target, 

Apv = (19) 



where ct/j^l = {d?a / dVLdE')R^L. Since the purely electromagnetic and purely weak contribu- 
tions to the cross section are independent of electron helicity for <^ M|, they cancel in 
the numerator, leaving only the interference term. The denominator, on the other hand, 



contains all contributions, but is dominated by the purely e. 
terms of structure functions, the asymmetry can be written 



ectromagnetic component. In 
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A 



PV 



/ GfQ^ 



a 



■7Z 



(20) 



where g'\ = — 1/2 and Qy = — 1/2 + 2 sin^ 6^ are the axial and vector couplings of the Z 
boson to the electron, with 6w the weak mixing angle, and the functions Yi^s parametrize 
the dependence on y and on the W and ratios, 

l + {l-yY- y^[l + p2 _ 2pV(i?^^ + l)]/2 fl + R^^' 



>3 



1 (1 - - y^[l + p2 _ 2p2/(i?7 + l)]/2 
1 + (1 - y)2 - y^[l + p2 _ 2p2/(i?7 + l)]/2 \1 + Ry 



1 + 



(21a) 
(21b) 



In the limit of — t- 00, where p — )■ 1 and R^''^^ — )■ 0, the kinematical factors simplify to 
Fi ^ 1 and F3 ^ [1 - (1 - yf]/[i + (1 - y?]. 



1. Proton asymmetry 

The proton PVDIS asymmetry is shown in Fig. [T2] for Q"^ = 2 and 10 GeV^ in the form of 
the ratio of the target mass corrected to uncorrected asymmetries. For all prescriptions the 
TMC effects are maximal at x ~ 0.7, where they are of the order of 3 — 4% at = 2 GeV^ 
and ^ 1% at = 10 GeV^. The results are slightly higher for the ^-S and AQ corrections 
(which are virtually indistinguishable) than for the OPE and EFP (which are also almost 
identical). The small size of the effects is principally due to the strong cancellation of the 
TMCs in the Fi structure functions, namely, (F^^/F^)™^ ^ {F^^/F^Y^\ even though 
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FIG. 12: Ratio of target mass corrected (^pv) to uncorrected (^py) PVDIS asymmetries for the 
proton at Q2 = 2 GeV^ (left) and = iq QeV^ (right), for the OPE (sohd, red), EFF (short- 
dashed, blue), ^-S (long-dashed, green), and AQ (dot-dashed, orange) TMC prescriptions. Note 
that the AQ and ^-S results are almost indistinguishable, as are the EFP and OPE prescriptions. 

|^TMC^^^(o)| ^ I gj^ high X. Overall, the results indicate that the asymmetries themselves 
are less sensitive to TMCs than are the LT ratios W'"''^ on which the asymmetries depend. 

Since one of the main goals of the proton PVDIS measurements will be to reduce the 
uncertainties on PDFs at large x, particularly on the d/u ratio, it is instructive to com- 
pare the magnitude of the TMC effects with the expected sensitivity of the asymmetry to 
different possible PDF behaviors at large x. In Fig. [13] we show the proton asymmetry 
y4pv computed from the full range of CJ PDFs [3] including minimal and maximal nuclear 
corrections (shaded bands) relative to the central PDF fits. The uncertainty band increases 
with increasing x, reflecting the larger uncertainty on the d quark PDF at large x, and in 
the absence of TMCs ranges from 3% at x = 0.6 to ~ 11% at x = 0.8 for both = 2 
and 10 GeV^. This is signiflcantly larger than the TMC uncertainty band in Fig. [T21 where 
the spread of the TMC model predictions is ^ 1%, even though the absolute target mass 
effect is somewhat larger. 

The effect of TMCs on the PDF uncertainty, illustrated in Fig. [13] for the OPE prescrip- 
tion, is to reduce the uncertainty band at large x for the lower value, in analogy with 
the effect seen in Fig. [7] for the Rnp ratio, with strength moving from lower x to higher x by 
the X — )■ ^ rescaling of the structure functions. The slightly different effects of TMCs on the 
various structure functions present in the asymmetry render the uncertainty band on Apv 
more asymmetric at = 2 GeV^. At the higher = 10 GeV^ value, the impact of TMCs 
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FIG. 13: Proton PVDIS asymmetry Apy at 



^2 GeV^ (left) and = 10 GeV^ (right) for CJ 

,(0) 



PDFs with minimal and maximal nuclear effects 71] (shaded bands) , relative to the asymmetry Apy 
for the central CJ PDF fits, with no TMCs (dashed, blue) and using the OPE TMC prescription 
(solid, red). 

on the uncertainty band is reduced considerably, with the two bands (corresponding to no 
TMCs and the OPE TMC prescription) approximately coinciding for all x. 

The conclusion from the combined results of Figs. [T2]and[T3]is that the effect of TMCs and 
particularly their uncertainties can be minimized in the Apy ratio by measuring the asym- 
metry at values of ~ 10 GeV^ or higher; at lower Q^, although the TMC uncertainties 
are not large, some residual corrections will need to be applied in the range 0.4 ^ x ^ 0.9, 
where the TMCs are ~ 1% or higher. 



2. Deuteron asymmetry 



Unlike for a proton target, for PVDIS on an isoscalar deuterium nucleus most of the 



dependence on PDFs cancels i: 



one assumes that PDFs in the proton and neutron are 



related by charge symmetry 25|]. In fact, in the valence quark region {x > 0.5) where sea 



quarks and gluons can be neglected, the deuteron asymmetry can be written at leading order 



as 
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30 



Ad ^ 



I ^^^^^^^ ~ + ys9v{29l - 9'a)) , [x » 0], 



(22) 



where = -1/2 + (4/3) sin^ ^^y, 9v = 1/2 - (2/3) sin^ ^^y, 9l = 1/2, and = -1/2. 
(Note that the conventions for gy^^ in Ref. |23| differ by a factor 2 compared with those used 
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here, although the asymmetry is of course independent of the convention.) Consequently 
accurate measurement of deuteron PVDIS has been proposed as a sensitive test of either the 
weak mixing angle sin^ Oy/ (deviations of which from its Standard Model value may signal 
the presence of new physics), or more conventionally of charge symmetry violation (CSV) 
in PDFs. 

Nonzero values of 5u and 5d are predicted in nonperturbative models of the nucleon to 



arise from quark mass differences and electromagnetic effects (for a review see Ref. 46|). 



and can also be generated from radiative QED corrections in evolution 49l-l5l|. Defining 
charge symmetry violating PDFs by 

5u = uP- cT, 5d = (F- m", (23) 

the PVDIS asymmetry (122!) in the presence of CSV is modified according to 

(2^7y,A - 9va) ^ i'^dvA - 9vA)i'^ + Av,a), (24) 
where the fractional CSV corrections are given by 



These approximate expressions serve to illustrate explicitly the role of CSV in the PVDIS 
asymmetry; in practice, however, the full deuteron asymmetry can be computed including 
the effects of CSV at NLO, as well as sea quarks and gluons. 



Using the MRSTQED parametrization of PDFs [SOj, which generates nonzero values of 
6u and 6d through radiative QED effects, the effect of CSV on the deuteron asymmetry Apy 
is illustrated in Fig. [TH In the valence quark region the CSV effect is small at intermediate 
x, x ~ 0.4, but increases to around 1% at x ~ 0.8. This is roughly comparable to the earlier 



fit in Ref. j49| which parametrized the CSV PDFs as 6u - 6d = 2k^{\ - xf{x - 0.0909), 
with K = —0.2 as the best fit parameter. (The constraints on k, were found to be relatively 
weak, however, and values of k, from —0.8 to +0.65 produced fits at the 90% confidence 
level, with effects on the asymmetry ranging from ^ 4% to 8% over the range 0.4 ^ x ^ 
0.8.) Deviations of the full NLO result from the valence approximation appear already 
at X ^ 0.7, however, and these differ quite markedly at small x, as Fig. [H] indicates. 
Interestingly, the full asymmetry becomes larger at smaller x because of CSV effects in 
the light sea quarks, which produce an asymmetry of about 2% at x ~ 0.2. On the other 
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FIG. 14: Deuteron PVDIS asymmetry including CSV effects, relative to the asymmetry with 
charge-symmetric PDFs, at = 2 GeV^ (left) and = 10 GeV^ (right). The CSV PDFs are 



computed from the MRSTQED parametrization 



50( 1 for the full asymmetry (solid, red) and for the 



Ml 



valence approximation (double-dashed, black), and from the K-dependent fit (see text) in Ref. 
(long-dashed, green). The effects of TMCs on the full asymmetry with the MRSTQED PDFs are 
illustrated for the OPE prescription (short-dashed, blue). 

hand, cleanly separating the CSV effects from sea quark and gluon contributions, which 
introduce additional x dependence beyond that in Eqs. fl22l) . fl24|) and fl25|) . as well as 
possible differences between CSV in valence and sea quark PDFs, becomes more challenging 
at small x. 

With sought-after CSV effects that could be ^ 1 — 2%, it is vital to quantify the impact 
of TMCs on the deuteron PVDIS asymmetries and in particular the TMC prescription 
dependence. The effect of TMCs on the full asymmetry relative to the charge-symmetric 
asymmetry is negligible at a; ^ 0.5, but decreases the CSV signal by up to 50% at a; ~ 0.8, 
as Fig. [m demonstrates for the OPE prescription. The model dependence of TMCs is 
illustrated for the various prescriptions in Fig. [T5|, where the ratio of asymmetries is shown 
with TMCs to those without TMCs. The net effect is very small, peaking at ~ 0.1% at 
X ~ 0.4, even at the = 2 GeV^ value. The TMC prescription dependence of this ratio 
is even smaller, making it essentially negligible on the scale of a CSV signal of ~ 1%. If 
the target mass corrected asymmetries were calculated with the charge symmetry violating 
MRSTQED PDFs, the effect would be somewhat larger, peaking at ~ 0.3% around x ~ 0.4. 
However, the TMC model dependence is still negligible at around 0.05%. As expected, 
the impact of TMCs on the deuteron asymmetries at the larger = 10 GeV^ value is 
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FIG. 15: Ratio of target mass corrected to uncorrected PVDIS deuteron asymmetries Apy at 
Q2 = 2 GeV2 (left) and = iq GeV^ (right), for the OPE (solid, red), EFP (short-dashed, 
blue), ^-S (long-dashed, green), and AQ (dot-dashed, orange) TMC prescriptions. 



considerably smaller. It is therefore likely that TMCs would only play a role in deuteron 
PVDIS measurements if the CSV effects were on the scale of a fraction of a percent, at 
which point they would not be discernible within the expected precision of the experiment 
[isl . The corrections due to nuclear smearing in the deuteron would have similarly negligible 
effects on the measured deuteron asymmetry. 



IV. CONCLUSIONS 



With the increased precision and kinematic reach of new experiments planned in the 
next few years, particularly at Jefferson Lab with its 12 GeV upgrade, the need for reliable 
theoretical tools with which to analyze the data is becoming ever more pertinent. This is 
especially true for data that will be taken at large values of a;, where a number of different 
subleading effects come to the fore. In this work we have performed a comprehensive analysis 
of one class of such corrections, namely those associated with finite values of x^M^/Q^, 
or target mass corrections. We have detailed several approaches to TMCs, including the 
standard OPE method, as well as prescriptions based on collinear factorization, and have 
compared their effects on various spin-averaged structure functions at next-to-leading order. 

For the TMCs computed via the OPE, we find that the l/Q^ and approximations 
to the full results are accurate only up to a; ~ 0.6, beyond which the series displays rather 
slow convergence. Such an expansion has been proposed to avoid the threshold problem at 
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X — 1; our findings suggest, however, that a low order expansion may not be apphcable as 
X — > 1. Numerically, we find that TMCs in the OPE approach are very similar to those 
computed via the EFP implementation of coUinear factorization, especially for the vector 
structure functions Fi (or F^) and F2. This can be demonstrated analytically, through the 
equality to order of the prefactors associated with the leading terms. The comparative 
phenomenology of these prescriptions has not previously been addressed in the literature. 

Similarly, the ^-scaling and AQ prescriptions, which are derived from different approx- 
imations within the coUinear factorization framework, yield corrected structure functions 
that closely track each other over much of the x range accessible experimentally. In all cases 
the magnitude of the TMCs, and in particular their model dependence is, not surprisingly, 
significantly more important at low values (~ 2 GcV^). Target mass corrections are 
suppressed with increasing Q^, although even at ~ 10 GeV^ they are not negligible for 
some observables at very large x. The greatest model dependence of TMCs arises for the 
longitudinal structure function, where because of the mixing between the Fl and F2 struc- 
ture functions the effects for the OPE and EFP prescriptions are significantly larger than 
for the AQ and ^-scaling approaches, where no mixing occurs. 

In addition to quantifying the impact of TMCs on structure functions, wc further dis- 
cussed the limitations these place on unambiguously extracting information on PDFs (such 
as the d/u ratio or charge symmetry violation) from observables. For the ratio Rnp of neu- 
tron to proton F2 structure functions we make the interesting observation that at low 
not only is one subject to greater TMC uncertainties than at large Q^, but the x — > ^ 
rescaling due to TMCs effectively also decreases the sensitivity to the d/u ratio at large x 
that measurements of Rnp attempt to constrain. 

For parity-violating DIS from the proton, the effects of TMCs and perturbative NLO 
radiative corrections are similar in both the electromagnetic and interference LT ratios 

and R'^^ , with < 4 — 5% differences for Q'^ — 2 GeV^ at intermediate and larger x. For 
the deuteron the differences between R'^ and RP'^ are smaller in the valence quark dominated 
region, with negligible dependence on the TMC prescription, but become larger at very small 
X (^ 8% and 4% at = 2 and 10 GcV^, respectively) through gluonic contributions at 
NLO. The magnitude of TMCs in the R'^^ ratio itself, however, is significant at large x, 
especially for the OPE and EFP prescriptions. We also considered the effects of nuclear 
corrections in the deuteron on the LT ratio, which become important for x ^ 0.4 — 0.5; 
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however, the similarity of these with the effects on the electromagnetic LT ratio leads to 
nuclear corrections largely canceling in the PVDIS asymmetry. 

The effects of TMCs on the parity-violating asymmetries themselves are generally rather 
small, especially at higher values, ~ 10 GeV^, although at lower some residual 
TMC dependence is evident in the case of the proton asymmetry. Measurements of the 
proton PVDIS asymmetry are planned to provide a unique combination of PDFs in order to 
constrain the dju ratio at large x 13|. For the deuteron, the size of TMCs is about an order 
of magnitude smaller than the expected CSV effects in PDFs, which are estimated to be at 
the 0(1%) level. On the other hand, while the corrections to the LT ratios and asymmetries 
computed here have been perturbative, nonperturbative effects such as those associated 
with nonzero parton transverse momentum in the nucleon can produce additional strength 
in the longitudinal structure functions 3]. This may be particularly relevant for the ratio 
W'^ ^ whose phenomenology is essentially unknown at low . Estimates of nonperturbative 
contributions to E?^ would therefore be necessary before making more definitive conclusions 
about its role in PVDIS. 

In the future, additional effects not discussed here may need to be considered at large 
X, principal among which are dynamical higher twist corrections associated with nonpertur- 
bative multi-parton correlations. These are very difficult to compute from first principles, 
and only rudimentary model estimates have been available to date. Further insight into the 
relation between TMCs and higher twists may also shed light on the threshold problem, 
whereby the target mass corrected structure functions remain finite at the nucleon elastic 
scattering point, X — X, cis well as on the difference between the various TMC prescriptions. 
Other corrections that may affect future analysis of large-x data are threshold resumma- 
tions, which involve formally summing, to all orders in a<(, terms containing logarithms of 
1 — X that become large 1. The results on the phenomenology of the target mass 

corrections contained in the present work should provide a benchmark for future theoretical 
and experimental investigations of these additional corrections. This analysis can also be 



extended to the spin-dependent sector [53 



56l ]. where the phenomenology of the coUinear 



factorization framework in particular has not been as fully developed. 
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